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Abstract. We investigate to what extent 18 Ne can be descibed as a three-body system made of an inert 
le O-core and two protons. We compare to experimental data and occasionally to shell model results. 
We obtain three-body wave functions with the hyperspherical adiabatic expansion method. We study the 
spectrum of 18 Ne, the structure of the different states and the predominant transition strengths. Two + , 
two 2 + , and one 4 + bound states are found where they are all known experimentally. Also one 3 + close 
to threshold is found and several negative parity states, 1~ , 3~, 0~, 2~, most of them bound with respect 
to the ls O excited 3~ state. The structures are extracted as partial wave components, as spatial sizes of 
matter and charge, and as probability distributions. Electromagnetic decay rates are calculated for these 
states. The dominating decay mode for the bound states is E2 and occasionally also Ml. 
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1 Introduction 

Nuclear cluster structures appear in various disguises es- 
pecially in light nuclei. The cluster constituents are often 
nucleons and cv-particles, possibly combined with a core- 
nucleus. These structures, which appear both as ground 
and excited (resonance) states, are sometimes well de- 
scribed as three-body systems. The conventional wisdom 
is that prominent clusters are most likely to appear close 
to the threshold energy for fragmentation into the cluster 
constituents. This implies that cluster structures for or- 
dinary bound nuclei are more likely to appear in excited 
states than in ground states, except for dripline nuclei 
where the ground state is close to the nucleon threshold 
and dominating one or two-nucleon structures appear [TJ 

Well-known three-body examples are the first + reso- 
nance in 12 C [5], the lowest + and 2 + states in 6 He, 6 Be 
and 6 Li [B], the 3/2~ ground state in n Li [7], three excited 
bound states in 12 Be [S], the ground state and four reso- 
nances in 17 Ne [3], the ground state and several resonances 
in 9 Be [10] . and three resonances in 5 H [IT] . Other nuclear 
states have significant admixtures of non-cluster structure 
( 12 C(2 ± )) [I2J while some states are far better described 
without any cluster structure ( 12 C(1 ± )) [5]. One line of 
investigation is to carry out the three-body computation 
for a given system and compare the computed observables 
with known measured values and then predict others. 

If the computed bulk structure of a nuclear state matches 
measurements the description is an immediate success. 



However, even for cases where no traces of any three-body 
structure can be found the computation can be consid- 
ered a necessary ingredient to describe the three-body de- 
cay of an underlying many-body resonance, examples are 
12 C(1 ± ) [13]. Resonances decaying into three clusters are 
now investigated accurately in details in complete kine- 
matics [14] . Both structure and dynamic evolution from 
small to large distances are important in a description of 
the momentum distributions of the decay products. A few 
such decaying structures have been investigated theoreti- 
cally and compared to available data jTSJUD] ■ 

Recently the decay products, two protons and 16 0, 
from the 6.15 MeV 1~ state in 18 Ne was measured [171116] . 
In the same nucleus the 4.522 MeV 1~ state received at- 
tention as a doorway state to produce the water molecule 
18J which has the same number of neutrons, protons and 
electrons as the 18 Ne-atom. Also the first 3 + state of this 
nucleus plays an important role in astrophysics related 
to the abundance ratio between 18 F and 17 F [19]. The 
17 F(p, 7) 18 Ne cross section has been investigated in a two- 
cluster coupled-channel model [20] , and recently also mea- 
sured directly 21 . The first question in this connection is 
obviously which structures have these, and perhaps other, 
18 Ne states. The low-energy spectrum of 18 Ne is typical 
for a quadrupole vibration with an equidistant spacing be- 
tween 0+, 2+, and a triplet of 0+,2+,4+ states [22]. Still 
higher at and above the two-proton threshold a + and 
a 3 + appear with a number of other states without an 
obvious recognizable pattern, see FigJT] 
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Fig. 1. Measured spectra compared to computations with the 
interactions I,II,III (see text and table [1} relating to the + 
ground state of the core and IV relates to the 3~ core exci- 
tation. The threshold of E = corresponds to infinite zero 
separation of the two protons and the oxygen core. 



The different states may originate from separate struc- 
tures described for example as vibrations, single-particle 
excitations, pairing correlations, or obtained in combina- 
tions of even more complicated few or many-body features. 
In particular the coupled two-body cluster model provides 
one type of structure information |20j . Also a number of 
interacting shell model calculations have provided struc- 
ture information about the low-lying excited states of the 
A = 18 isobaric system, 18 18 F 18 Ne, see e.g. [25P3]. 
The results are in general that many of the states are 
more complicated than two nucleons and the 16 0-core. 
This is not surprising when the excitation energy is suf- 
ficiently large to accomodate core excitations. However, 
the shell model is designed to describe spatially confined 
bound state structures without strong cluster configura- 
tions beyond the chosen core- valence division. This means 
that large-distance structures are inaccessible or inaccu- 
rate in shell model calculations. This applies in particular 
to resonances and doorway states in reactions. 

Specifically, three-body decays of resonances cannot be 
described by two-body cluster or shell models. The reso- 
nance structures necessarily change from the many-body 
short-distance behavior to three-body clusters at interme- 
diate distances. In addition, the two-body structure is in- 
adequate and the three-body structure itself often change 
dramatically from intermediate to large distances [FjllfjlfTUl 
[T2] . The three-body structures must be accurately de- 



scribed to meet requirements of up-to-date measurements. 
In other words shell model calculations necessarily must 
be supplemented by few-body calculations as provided in 
the present paper. 

To clear the road towards computing the three-body 
decays measured in [T71ITB] we start with assuming few- 
body structures to see how far this will bring us in a 
quantitative description of the various states. Since 17 Ne 
( 15 0+p+p) is Borromean an extra neutron suggest a four- 
body structure but a neutron and 15 form a strongly 
bound doubly magic nucleus, 16 0, and a three-body struc- 
ture of 16 0+p+p is probably a better starting point. Un- 
fortunately it is then unlikely that the 1~ states simulta- 
neously are simple three-body structures maintaining an 
16 ground state core. At least the 3~ state in the 16 
core [22] can be expected to contribute. 

In the present paper we attempt to describe the low- 
lying (bound and resonance) states in 18 Ne as three-body 
states. If possible this is a huge simplification from the full 
problem of 18 interacting nucleons. These investigations 
are a generalization of the classical nucleon-core model 
and its extension to two mutually interacting particles oc- 
cupying single-particle levels provided by a core but with- 
out additional nucleon-core interaction. In any case, the 
results are a prerequisite for description of three-body de- 
caying resonances like the measured 1~ state [IrJlll7| . The 
paper is organized as follows. In section [2] we briefly de- 
scribe the notation by sketching the three-body method 
and the constraints used to determine the crucial two- 
body interactions. The structures are shown in section [3j 
and the sizes and electromagnetic transitions are given in 
|U Finally, section [5] contains a summary and the conclu- 



2 Method and interactions 

The principal model assumption is that 18 Ne can be de- 
scribed as a three-body system made by a 16 core and 
two protons. The wave functions for the different bound 
states are obtained with the hyperspherical adiabatic ex- 
pansion method. A detailed description of the method can 
be found in l25l. 



2.1 Theoretical formulation 

This method solves the Faddeev equations in coordinate 
space. The wave functions are computed as a sum of three 
Faddeev components tp^(xi, yj (i=l,2,3), each of them 
expressed in one of the three possible sets of Jacobi co- 
ordinates {a?i,yj}. Each component is then expanded in 

terms of a complete set of angular functions {<fin^} 

^ = 4j2 E/^n^A);^ = {a h a Xi ,n yi }), 

^ n 

(i) 

where p = \J x 2 + y 2 , cti — arctan^/j/i), f2 Xi , and Q Vi are 
the angles defining the directions of Xi and y^. Writing the 
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Faddeev equations in terms of these coordinates 
be separated into angular and radial parts: 
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where Vjk 
and k, A 2 



is the two-body interaction between particles j 
is the hyperangular operator [25 and m is the 
normalization mass. In Eq.@ E is the three-body energy, 
and the coupling functions P nn > and Q nn > are given for 
instance in [35] ■ The potential Vzb{p) is used for fine tuning 
to take into account all those effects that go beyond the 
two-body interactions. 

It is important to note that the set of angular functions 
used in the expansion ([I]) are precisely the eigenfunctions 
of the angular part of the Faddeev equations. Each of them 
are in practice obtained by expansion in terms of the hy- 
perspherical harmonics. Obviously this infinite expansion 
has to be cut off at some point, maintaining only the most 
essential components. Specifically the contributing partial 
waves increase with energy and distance. We include suf- 
ficiently many higher partial waves to reach convergence. 

The eigenvalues X n (p) in Eq.@ enter in the radial 
equations ([3]) as a basic ingredient in the effective radial 
potentials. Accurate calculation of the A-eigenvalues re- 
quires, for each particular component, a sufficiently large 
number of hyperspherical harmonics. In other words, the 
maximum value of the hypermomentum (K max ) for each 
component must be large enough to assure convergence of 
the A-functions in the region of p- values where the f n {p) 
wave functions are not negligible. 

Finally, the last convergence to take into account is the 
one corresponding to the expansion in Eq.([T]). Typically, 
for bound states, this expansion converges rather fast, and 
usually three or four adiabatic terms are already sufficient. 



where i is the relative orbital angular momentum between 
the two protons, and s = Si + s 2 is the total spin. The 
strengths are in MeV and the ranges in fm. We shall refer 
to this potential as the gaussian proton-proton potential. 
In actual three-body computations we have tested, see e.g. 
[HIE]) by using other nucleon-nucleon potentials like the 
Argonne and the Gogny potentials. The three-body results 
were always indistinguishable. 



2.3 Proton- 16 potential 

The other two-body interaction is related to the proton- 
16 system. The core, 16 0, has intrinsic spin and parity 
+ and the proton has spin 1/2. The most general spin 
dependence is then of spin-orbit form and each orbital 
angular momentum potential has the form 



V (e Hr) = V®{r) + V®(r)l ■ s p + V c (r) 



(5) 



where t is the proton-core relative orbital angular mo- 
mentum and s p is the spin of the proton. These potentials 
should be parametrized to reproduce low-energy scatter- 
ing properties. We assume gaussians for the radial shapes 



of all terms, i.e. central, Vc (r), and spin-orbit, V}o'(r). 
The range b of the gaussians has to be related to the size 
of the 16 0-core. We choose b = 2.60 fm, as selected in 9 
for the proton- 15 potentials. The strengths are then left 
as adjustable parameters. The Coulomb potential, Vc(r), 
can be obtained either from 16 O as a point particle or with 
a gaussian charge distribution. For our purpose it suffices 
to use the potential from a point charge as we did for the 
proton-proton Coulomb interaction. 

The most prominent features in low-energy scattering 
data are reflected in properties of bound states and reso- 
nances, where the dominant features in turn are energies 
of these structures. We shall therefore first aim at repro- 
ducing these energies. The two-body system is 17 F which 
has two proton bound states, i.e. a l/2 + -state at —0.105 
MeV and a 5/2 + -state at —0.600 MeV measured relative 
to the two-body threshold [22 . The s-wave has no spin- 
orbit term and the strength of Vc °^ (r) can be determined 
to reproduce the energy —0.105 MeV. This should be the 
potential reproducing the experimental s- and p-wave proton-second s-state as the first is occupied by core protons and 



2.2 Proton-proton interactions 

The two-body low-energy scattering properties are cru- 
cial in the description of weakly bound systems. The de- 
tailed short-distance behavior is relatively unimportant. 
Thus we adjust the parametrized two-body interactions 
to known low-energy properties. In the present case this 
means the nucleon-nucleon interaction or, to be specific, 
the proton-proton interaction. We use a simple short-range 



proton scattering lengths and effective ranges. This as- 
sumes that effects of the Coulomb interaction are removed 
from these scattering parameters. Obviously the Coulomb 
potential is then afterwards added in the final potential. 
We assume the protons are point-like particles and the 
Coulomb potential is then e 2 /r. 

The short-range nucleon-nucleon potential contains cen- 
tral, spin-orbit (£ ■ s), tensor (512) and spin-spin (si ■ S2) 



consequently Pauli forbidden. To exclude the lowest s- 
state in three-body computations we can either use a shal- 
low potential with only one bound s-state at —0.105 MeV, 
or construct a phase equivalent potential (P.E.P) with one 
less bound s-state [27] ■ For the d-state we keep the same 
central potential strength we use for s-states while adjust- 
ing the strength of the spin-orbit term to give a d 5 /2-statc 
at -0.600 MeV. 
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Fig. 2. Differential Cross Sections for 16 0(p,p) 16 0. The com- 
puted lines were calculated using s and d-waves generated with 
the potentials in table [T] 



These potentials now also lead to elastic cross sections, 
or equivalently, phase shifts for each set of quantum num- 
bers. We compare in table [2] the computed phase shifts 
with measured values from [35]. The Si/ 2 and d 5 / 2 phase 
shifts are matching the data perfectly as expected because 
the positions of the bound states are well determined in 
our fits to match the measured values. On the other hand 
the d 3 / 2 phase shifts deviate by several degrees although 
both calculated and measured values are small. This par- 
tial wave is unimportant for the low-energy structures, 
and we have not attempted any adjustment to these ob- 
servables. We also computed the differential cross section 
as measured for several angles in [29,30,31,32 . As we can 
see in Fig. [2l they are in perfect agreement with results 
from calculations with our potentials including only s and 
d-waves for energies up to 1 MeV where the higher par- 
tial waves begin to contribute. This is sufficient as the 
present computations almost exclusively only need s,p 
and d- waves. If more waves occasinally are needed we use 
the same potential parameters as for the d-waves. 

To determine the p-wave two-body interaction the usual 
procedure would be to reproduce negative parity 1/2 or 
3/2-states. Such two states are found above threshold in 
17 F at 2.504 MeV and 4.04 MeV for 1/2 or 3/2, respec- 
tively, see j33ll34ll35] . However, the sequence is opposite 
the established order from the spin-orbit splitting. Fur- 
thermore, these p-states should then correspond to single- 
particle excitations into the p — f shell which first should 
appear at substantially higher energies. Two choices seem 
at first to be possible for the p-wave interaction. The first 
is to enforce a p- wave potential to reproduce api/ 2 -energy 
at 2.504 MeV with the opposite sign of the spin-orbit po- 
tential perhaps with a strength related to the p3/ 2 -energy 
at 4.04 MeV. The second is to believe that these observed 
negative parity states are complicated many-body states 
without any influence on the three-body structure. This 
could be implemented by using the established shallow 
s-potential with the spin-orbit term from the d-wave po- 
tential. 

There is also a third option, which relates these levels 
to the 3~ excited state of the 16 0-core at 6.13 MeV above 



Table 1. Strengths in MeV of the central, V} '(r) = 

^e-W^'', spin-orbit, V s ( P(r) = S { a l J e^^, and spin- 

spin, Vaa ( r ) = Ssi* 1 e~ < - r ^ 1 ' 33 ' , potentials entering in Eq.© for 
the four interactions used in the calculations. The ranges are 
b c = b so = b ss — 2.60 fm in all cases. For the deep s-wave 
potentials the lowest s state is removed by construction of the 
phase equivalent potential (P.E.P.) from the central strength 
of the d-wave potential. The last column refers to potentials 
built on the 3~ excited state of 16 O. 
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Table 2. Phase Shifts for s and d- waves of F. For each energy 
the first row shows the experimental values from [2S] and the 
second row gives the computed phase shifts with potential I. 
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the ground state. This core-state coupled to a 5/2 + proton 
single-particle state could produce the two observed l/2~ 
and 3/2~ states. A similar coupling to a l/2 + proton state 
produce 5/2~ and 7/2~ states. The core excited 1 _ state 
at 7.12 MeV could also couple to the 5/2 + and l/2 + single- 
particle states to give the 3/2~, 5/2~, 7/2~ and the l/2~, 
3/2~ states. Then the lowest-lying l/2~ state would not 
be related to the lowest d-wave but to the higher lying s- 
wave. The 3/2~ would be a mixture of both s and d-waves. 
We restrict ourselves to explore the simplest combination 
relating to the 3~ excited core state. 

These negative parity states could then have contri- 
butions from both the core excited 3~ state and the + 
ground state of the 16 core. To the degree that they are 
decoupled in the two-body states they would also be de- 
coupled in the three-body states. Furthermore, the lowest 
possible partial waves for most of the low-energy three- 
body J w states correspond uniquely to either 3 _ or + 
core states. Other contributions are less favored by either 
relative energy or core excitation. Thus decoupling on the 
three-body level could be rather well fulfilled. 

We have now established several options for the two- 
body proton-core potential. The strengths of the resulting 
different potentials are given in table[TJ We include specific 
interactions for s, p, and d-waves. In all cases are the two 
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bound state energies reproduced. Potentials I and II em- 
ploy the same deep potential for both s and d- waves, the 
phase equivalence is used for s- waves. Potential III main- 
tain the d-wave from I and II, whereas a shallow potential 
with one bound state is used for s-waves. The p-wave po- 
tentials I and III use the central s-potential from III and 
the spin-orbit from I and II. In II the p-wave is adjusted 
to give the measured pi/ 2 energy at 2.504 MeV. The most 
likely candidate as spin-orbit partner of the 1/2™ state is 
a 3/2" resonance at 4.04 MeV (above threshold). There- 
fore we adjust the central and spin-orbit depth to fit these 
energies while keeping the same configuration as I for s 
and d- waves. 

Potentials I, II, and III are based on the spin zero 
ground state of 16 0. Potential IV is based on the 3~ ex- 
cited state of 16 O where the d-wave is the decisive compo- 
nent in the description of the low-lying states. This finite 
spin of the core requires the generalization of the potential 
to include the spin-spin term, i.e. 

V^(r) = V^\r) + V<?(r)t ■ s p + v£{r)a c ■ j p + V c (r), 

(6) 

where we maintain the gaussian shapes of range 2.60 fm 
of the radial potentials, and the Coulomb potential again 
is for point particles as for potentials I, II and III. We use 
the energies of the 1/2™ and 3/2™ states to determine the 
central and spin-spin strengths for the ds/2-wave coupled 
to 3~ of the core. This leaves the spin-orbit strength as 
a free parameter which is chosen to be similar in size to 
the values of the other potentials. It only has to be large 
enough to place the d 3 / 2 -waves at sufficiently high energies 
to make their effects negligibly small. 

Coupling of Si/2-waves and 3™ produce 5/2™ and 7/2™ 
states which however also arise from the d 5/ / 2 couplings. 
These relatively high-lying states are expected to con- 
tribute very little to the low-lying three-body states. We 
use the energy, 3.257 MeV above threshold, of the 5/2™ 
resonances to estimate the central strength of the s-wave 
potential while maintaining the spin-spin strength derived 
from the d-waves. Any p-wave would in principle mix into 
the positive parity states but both angular momentum and 
energy indicate negligibly small effect. We have kept the 
same potential as for s-waves adding the spin-orbit from 
d- waves. 

We have not attempted to reproduce the widths of 
these resonances as that would require more parameters 
like variation of the range of the potentials. 



3 Structure 

The angular eigenvalues obtained from eq.Q enter in the 
coupled set of radial equations eq. ([3j> as the crucial ingre- 
dient of the effective potentials V^J(p) — \ (A n (p) + -^). 
To solve the angular part of the Faddeev equations we use 
expansion of each Faddeev component on hyperharmonic 
wave functions. The angular wavefunctions are expressed 
in the Jacobi coordinates, x and y, where the correspond- 
ing orbital angular momentum quantum numbers, t x and 
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P (fm) p (fm) 

Fig. 3. The four lowest adiabatic effective potentials V^fl for 
positive parity states in 18 Ne built on the + core structure 
and with total angular momentum J n — + ,2 + ,3 + ,4 + and 
also ,V — 0™,3~ for 3™ core. The proton-core potential is I 
for the + core cases and potential IV for the 3™ core case as 
given in table [T] 

£ y , couple to the total orbital angular momentum L. The 
parity 7r is then given by the odd or even character of 

i x + £ y . The spins of the two particles connected by the 

x coordinate couple to s x , that in turn couples with the 
spin of the third particle to the total spin S. Finally L and 
S couple to the total angular momentum J of the three- 
body system. The last quantum number of the basis is 
the hypermomentum K — In + i x + l y where the non- 
negative integer n is the number of nodes. For each set of 
angular quantum numbers we include all K- values smaller 
than a given K max chosen to guarantee convergence for all 
necessary hyperradii. 



3.1 Effective potentials 

The potentials in eq.Q determine the structure of the 
18 Ne states. We show in Fig. [3] the lowest effective poten- 
tials for the bound states and resonances selected among 
the different structures we have investigated, i.e. J 7 ' — 
± , 1 ± , 2 ± , 3 d1 , 4 + built on both + and 3™ core-structures. 
We see attractive pockets in the two lowest potentials for 
the positive parity states built on the + core. These po- 
tentials are strongly influenced by the bound s and in- 
states in the proton-core potentials. At larger distances 
they asymptotically approach the energies, —0.6 Mev or 
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—0.1 Mev, of the 17 F bound states. This reflects the large- 
distance structure of one proton far away from the 17 F nu- 
cleus in the corresponding bound states, i.e. ground state, 
(^5/2, and first excited state, s\/2- 

For the total angular momentum, + , the last proton 
is in either d or s-waves around the corresponding 17 F 
structures. The 2 + potentials can have 17 F nucleus in the 
d 5 / 2 -state surrounded by the second proton with an an- 
gular momentum of either 0,2 or 4 in total allowing 5 
different couplings. We only include the three lowest of 
these potentials which approach —0.6 MeV. The fourth 
of the 2 + potentials in Fig. [3] approaches —0.1 MeV. It 
corresponds to the 17 F nucleus in the si/2-state where the 
second proton has angular momentum 2 around 17 F sub- 
system. At small distances two of these potentials exhibit 
rather appreciable attraction. 

The 4 + potentials allow proton angular momenta of 
2,4,6, and 4 around the 17 F nucleus in d or s-states, re- 
spectively. We find that two of these potentials approach 
—0.6 MeV and one approaches —0.1 MeV at large dis- 
tance. For the 3 + potentials the angualr momentum com- 
binations are 2,4,6, or 2,4, for the bound d and s-states, 
respectively. The lowest potentials again approach —0.6 MeV 
and —0.1 MeV. Most of the higher-lying potentials ap- 
proach zero reflecting a genuine three-body continuum 
structure. 

The lowest effective potentials for negative parity states 
built on the + core structure can also be computed with 
the help of the single-particle p- waves. Using p-wave inter- 
action from the s and d-waves in potential I we find again 
relatively attractive pockets in the lowest adiabatic poten- 
tials for 1~ and 3~ but they are almost totally absent for 
0~ and 2~. The large-distance approaches are found in 
all cases towards the s and d-wave two-body structures. 
Bound states or resonances of corresponding structures 
may then arise for 1~ and 3~. 

The potentials for the 0~ and 3~ states built on an 
excited 3~ 16 0-core state are also shown in FigEJ Simi- 
lar but less attractive potentials are found for 1~ and 2~. 
The threshold energy for all these states is then with re- 
spect to this core excited state at 6.13 MeV above the + 
ground state of 16 0. The two-body states used to adjust 
the interactions are bound with respect to this core excited 
state. These three-body potentials also exhibit attractive 
pockets at small distances. Their large-distance asymp- 
totics also reflect these two-body "bound states" where 
the three lowest potentials approach —3.6 MeV, —2.9 MeV 
and —2.1 MeV. These values correspond to the proton- 
core "bound" states of (£,S) = (2, 5/2), (0, 5/2), (2, 7/2) 
where I is the orbital angular momentum and S the total 
spin quantum number including the 3 _ from the core. 

At small distances several attractive pockets appear. 
For 0~ only one deep and broad potential can bind with 
respect to the 3~ excited core state. This is essentially due 
to an even combination of the (£, S) — (2, 5/2), (0, 5/2) 
partial waves. For 1~, 2~ and 3^ two attractive poten- 
tials are found with varying relative depths. They are 
mostly constructed from those combinations of (£, S) = 



Table 3. Spectrum of the positive parity bound states in 18 Ne 
for the different proton-core interactions in tableQ] The numer- 
ical results have been obtained without inclusion of a three- 
body potential in Eq.@. The two-proton separation energies 
are given in MeV. The last column gives the available experi- 
mental energies from |22| . 
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-1.38 
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-0.08 


-0.11 
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+0.04 



(2, 5/2), (0, 5/2), (2, 7/2) which allow spatially overlapping 
antisymmetric two-proton states. 



3.2 Three-body energies 

For each set of adiabatic potentials we solve the coupled 
set of radial equations in eq.©. The resulting eigenvalues 
are shown in tables [3] and E] for bound and unbound solu- 
tions, respectively. If the energies are decisive for applica- 
tions, as for breakup and decaying resonances, we can fine 
tune by use of the effective three-body potentials, Vsb(p), 
in eq.([3]). This would maintain the structure essentially 
completely unchanged. Such adjustments are not included 
in the eigenvalue tables. 

The unbound states are decaying resonances with a 
width arising from cluster, or equivalently two-proton, 
emission. To compute such continuum states complex scal- 
ing could be applied. States built on the excited 3~ core 
state with energies less than 6.13 MeV are bound states 
in the computation and complex scaling is not needed. 
They can only decay electromagnetically, or by the ne- 
glected coupling to the + ground state. In any case we 
focus in this paper only on the real part of the energies 
which above their respective thresholds are computed in 
two steps. First a sufficiently attractive three-body poten- 
tial is added to bind the state. Second the strength of that 
potential is varied and the resulting energy is extrapolated 
to the estimated energy obtained for zero strength. This 
is the so called analytic continuation of the coupling con- 
stant method ;36., which is especially relevant for states 
arising from the potentials for negative parity states based 
on the + core. 

The bound states shown in table [3] for different poten- 
tials are quite stable independent of choice of potential. 
In all cases the level ordering is reproduced and the ener- 
gies are also rather close to the measured values. We only 
adjusted the potentials to two-body bound state and reso- 
nance properties with the simplest possible radial shapes. 
Potential I leads in general to less binding but closer to 
measurements than the two other potentials. The devia- 
tions are less than 200 keV except for the 500 keV under- 
binding of the last + state. 

In potential I the phase equivalent s-wave potential 
produces a repulsive core at short distance and the va- 
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lence protons are pushed away from the center. This im- 
plies that potential III should bind more when proton-core 
s-wave configurations are substantial as for the first 2 + , 
the 3 + , and both the + states. The energies of both the 
second 2+ and the 4 + states do not depend on the cho- 
sen potential but matches almost perfectly the measured 
values. The 3 + state is weakly bound by less than about 
100 keV in the computation in agreement with the mea- 
sured value close to the breakup threshold. 

The results from potential I and II deviate surprisingly 
much for the two lowest states indicating that the p-wave 
components play a role. Then the description from po- 
tential II implies that these states have contributions of 
proton-core single- particle pi/ 2 character. This is inconsis- 
tent with the shell structure of the core-nucleus. We prefer 
potentials I and III with the much weaker p-wave attrac- 
tion. Potential I fits the experimental values better than 
potential II. 

As shown in table U we find a resonance at about 
3.10 MeV and a low-lying lj~ resonance at about 0.2 — 
0.4 MeV. Both states are based on the + core ground 
state. These values are more uncertain since they are ob- 
tained by extrapolation with a strongly attractive three- 
body potential. These structures are not present for 0~ 
and 2~ as already seen from the disappearance of attrac- 
tive pockets in the lowest adiabatic potentials for negative 
parity states based on the + core. 

An alternative to potential II in descriptions of the 
negative parity states is potential IV. We show several 
of the resulting energies in table |H In particular there 
appears a 1^~ state built on the excited 3~ core state and 
bound compared to this state by 3.98 MeV implying that 
it is an observable resonance at about 2.12 MeV. Since the 
energies of the two 1~ states differ by 2 MeV they may be 
decoupled in practice. The 0~ energy is about —5.45 MeV 
with respect to the core excited state and therefore at an 
energy of 0.58 MeV above the two-proton threshold. 

We find another 3~ "bound" state at about —3.09 MeV 
corresponding to 1.33 MeV. For 2~ we find two "bound 
states at about —2.92 MeV and —2.00 MeV corresponding 
to 3.41 MeV and 4.13 MeV. The second of these is a res- 
onance in the two-body continuum of the bound proton- 
16 0(3~) system and the other proton. There is no con- 
fining barrier and it easily leaks out corresponding to a 
large width. All these negative parity states can easily be 
matched to measured energies in the continuum. However, 
such a comparison is not very revealing due to the in- 
evitable inaccuracy from the three-body approximations. 
At least more structure information is needed. Still we 
show some of the lowest measured values in table |4] Many 
higher-lying levels are found experimentally. 



Table 4. Spectrum of the unbound, mostly negative parity, 
states in ls Ne for different proton-core interactions in table [T] 
The numerical results have been obtained without inclusion of 
a three-body potential in Eq.©. The two-proton separation 
energies are given in MeV with respect to the threshold for 
separation into two protons and the core in its ground state. 
These energies are given below the potentials used in the com- 
putation. The fifth and sixth columns give the available experi- 
mental energies and spin and parities, if known, from [25]. The 
computations with potential IV are for a core-excited state at 
6.13 MeV implying that these states with energies lower than 
6.13 MeV behave as bound states. 
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Table 5. Components included in the calculations for the 
states. The upper part correspond to the first Jacobi set (x 
between the two protons). The lower part corresponds to the 
second and third Jacobi sets (x from core to proton). The 6 th 
column gives the maximum value of the hypermomentum used 
for each component. The last three columns give the contribu- 
tion from each component to the + wave functions for poten- 
tials I, II and III respectively. Only contributions larger than 
0.1 is given. The two numbers for each component correspond 
to the Of and Ojj" states, respectively. 
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3.3 Wavefunctions 

The eigenfunctions are found as expansion coefficients on 
the hyperharmonic basis with quantum numbers for each 
of the Faddeev components, i.e. (£ x ,£ y , L, s x , S, J). Each 
cigcnfunction can be expressed in one set of Jacobi co- 
ordinates with corresponding probabilities depending on 



potential and quantum numbers. We shall in this section 
in details discuss the two bound + states and the 1~ res- 
onance located close to the threshold. For completeness 
we give the decompositions of the other states in the ap- 
pendix. 
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3.3.1 Bound states 

The available single-particle states for the two protons are 
d 5 /2 an d S1/2 orbits. Two-particle states of both protons in 
^5/2 produce the sequence of + , 2 + , and 4 + states where 
the odd angular momenta are forbidden due to the anti- 
symmetry requirement. Using two S\/ 2 states we can only 
produce a + state. One proton in each of the d 5 / 2 and 
S1/2 states produce one 2 + and one 3 + state. As seen in 
Fig.[T]they all appear in the computed spectrum. The par- 
tial wave decomposition reveal the microscopic structure 
of the states. 

We show the partial wave decomposition in table [5] for 
the two + states. The upper part using the first Jacobi 
system, where the ^-coordinate connects the two protons, 
exhibits only small variation between results from the dif- 
ferent potentials, the s-waves of about 80% dominate for 
both states in the three cases. 

The picture is very different in the proton-core Jacobi 
system where the variation with potential is larger. The 
2 state with potential I has about 74% in the Sy 2 config- 
uration which essentially means that both protons are in 
Si/2-states. With potential III the components d 5 / 2 and 
Si/2 are more even, i.e. about 55% d\, 2 in 0^ and 53% 
lS i/2 m ®2 ■ This reflects the lack of repulsion at short dis- 
tance in the s-wave interaction which favor s-waves in the 
ground state. With potential II we still get about 76% 
s\j 2 configuration in the 2 state. We also included com- 
ponents with i x and l y larger than 2 although their contri- 
butions are found to be negligible after the computations 
as we already mentioned in section f2. 31 

The configurations obtained here for the 2 state es- 
sentially only contains sd-waves. Early shell model cal- 
culations of both mirror nucleus 18 [35] and 18 Ne [33J 
gave about 33% and collective motion the remaining 67%. 
Our computed energy only deviates from measurements 
by about 0.4 MeV, see table[3J which is a typical deviation 
in such three-body calculations. This is therefore surpris- 
ing if 2/3 of the structure should have a completely dif- 
ferent origin. It is more likely that the sci-configurations 
contribute by more than 1/3 to the 2 structure. This 
may be reconciled with the shell model results if part of 
the collective motion also is of sd-character. 

The other four bound states of J 11 = 2 + ,3 + ,4 + are 
also decomposed in partial wave configurations and shown 
as tables in the appendix. Both 2 + and 3 + states consist 
of proton-core s and d-waves, and the 4 + state of solely 
d- waves. 



3.3.2 Unbound states 

With potential II it is a priori not excluded to find negative 
parity energies with resemblance to the measured spec- 
trum. The corresponding structures are on the other hand 
not expected to reproduce measured properties. The basic 
problem is the assumption of a single-particle p!/ 2 -state 
in the low-lying spectrum of 17 F. Potential IV, built on 



Table 6. The same as table [5] for the 1 1 state in 18 Ne with 
the + core. 
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as table [5] for the second 1 2 
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core excited state. 
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the core excited 3 _ state is in general expected to provide 
better structure properties. Then a number of 0~, 1~ , 2~, 
and 3~ states should arise as combinations of d and s 
states coupled to the 3~ core state. However, it is here 
worth emphasizing that there might be different, perhaps 
essentially uncoupled, structures of the same J T but built 
on different core states. We find such states with J*" = 1~ 
and 3~. 

The first l^f state is built on the + core ground state. 
The cluster model 14 0+a in [3D] cannot describe this state 
which is suggested as a candidate for burning 18 Ne into 
water while releasing a lot of energy [18| . The decompo- 
sition shown in table [5J is in the first Jacobi system seen 
to be dominated by s-wave components between the two 
protons. In the other Jacobi system both proton-core s 
and d-waves contribute about 20% and p-waves by twice 
that amount. In this way the attraction of the two interac- 
tions, S!/ 2 and g?5/ 2 , are optimized. The lowest adiabatic 
potential contributes by 93%. 

The second 1 2 state is built on the 3~ excited core 
state. The three-body decay of this state into two protons 
and the ground state of 16 O is measured and analysed 
in terms of sequential, virtual sequential and direct de- 
cay branching ratios [l"6lll~7] . Such decay must take place 
through couplings to other states as this state is bound 
with respect to the 3~ core excited state. The decompo- 
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sition in table [7J show dominance of p-waves in the first 
Jacobi system. In the other Jacobi systems this results in 
equal amounts of s and d-waves in the proton-core sub- 
system. This implies that the proton-core spin has to be 
5/2. The lowest and second potential contribute by 88% 
and 11%, respectively. 

We find several differences with respect to the shell 
model results in [24. where the structure of the different 
1~ states in 18 Ne are discussed in detail. The first five 
of these shell model 1~ states have for one choice of in- 
teractions respectively about 2%, 20%, 6%, 12%, 49% of 
configurations with 16 in the ground state. For compar- 
ison our two 1~ states contain 100% either ground or 3~ 
excited state of 16 O. The shell model results strongly in- 
dicate mixing of different excited states of 16 O coupled to 
the two protons. However, these shell model structures at 
most determine the short-distance behavior whereas the 
intermediate and large-distance structures can be com- 
pletely different and hence also the resulting momentum 
distributions after decay. 

The partial wave decomposition of all other computed 
unbound states are shown in the appendix. The 0~ and 
the two 2~ built on the core excited 3~ state all consist 
of proton-core s and d-waves. The two 3 _ states consist 
of proton-core p-waves and p and d-waves when built on 
the 3~ and + core states, respectively. 

The partial wave decomposition focus on the angular 
structure. The structures can be further characterized by 
the overlaps between valence wavefunctions of negative 
and positive parity states. Here it is necessary to remem- 
ber that the core structure differs, and true overlaps are 
zero. However, the valence part may have contributions 
of precisely the same partial waves which in turn has to 
be coupled to + or 3" to give the different total angular 
momenta. The overlaps can be estimated from the partial 
wave decompositions in the tables. These angular overlaps 
should be multiplied by radial overlap functions which in 
general are rather similar for low-lying states. The orbital 
and spin angular momentum couplings introduce in some 
cases another substantial reduction factor. It is in this way 
easily seen that the two + states and the 3~ state have 
the largest overlaps both exceeding 0.6. Also the 3 + and 
the 0~ states seem to overlap by more than 0.4 similar to 
the first 2 + and the first 2~ states. All other overlaps are 
rather small. 



4 Moments and transition probabilities 

The expectation values of the operators provide observ- 
ables for each state. The most interesting are those re- 
lated to sizes and lifetimes. We give details in the next 
two subsections. 



4.1 Relative sizes 

The sizes are observable quantities, where the simplest are 
the second moments of charge and matter distributions. 
These are the root mean square radii which often only are 



available for the ground state. The charge radius is most 
accurately obtained by electromagnetic probes like elec- 
trons. The matter radius is for light nuclei derived from 
measurements of interaction cross sections. The excited 
states are closer to the threshold for breakup and therefore 
more likely to develop a spatially extended halo structure. 
This would have observable implications for breakup cross 
sections. Effects of binding energy and angular momentum 
are both important [Tj . The three-body results are related 
to observables by including the finite extension of core 
distribution for matter and charge. In the present case we 
have for the matter distribution 

(r 2 ) = ^rl ore + ±(p*), (7) 

/ 2\ / 2 \ / 2 \ / 2 \ / 2 \ 

\P / ^vpp/ ~r -j^g Vc.pp) yj vpc) 13 vpjCp) i 

where r^ ore is the mean-square radius of the core. For the 
charge distribution we get 

/ 2 \ ^2 / 2 \ 136 , o ■ /o\ 

\T )ch — YQ T core,ch + ~^\ r pp) + 810 ,pp ' 

(r 2 )ch 9 ^lore.ch + <*jp)Ag + 17(7^)^/81 

(r 2 ) 10 r^ e + <^ p >/32+(r3 flp > oh /9 ' ( ' 

where r core ch is the mean-square radius of the core charge 
distributions. We assume that charge and matter distribu- 
tions are identical for the core. With these expressions we 
can always insert a different value for the core moments if 
better parameters become available. 

4.1.1 Bound states 

The root-mean-square, charge and matter, radii for the 
computed bound states are given in table [5] The results 
are essentially independent of the core-proton potential 
but the trends reflect that smaller binding energies give 
larger radii and viceversa. This is especially clearly seen 
for the J state where potential III gives more binding and 
smaller radius. Comparing the different + and 2 + states 
the tendency is also clearly that the smallest binding lead 
to the largest radius. The trends for angular momentum 
is that s-waves easier extend to larger distances whereas 
d-waves and higher are confined by the centrifugal barrier. 
This is seen for the 3 + state which is weaker bound with 
smaller angular momenta than the 4 + state, and conse- 
quently also significantly larger. 

When we assume equal matter and charge radius for 
the core we find that the charge radii are slightly larger 
than the corresponding matter radii. This is in spite of 
the "natural" reduction factor of 9/10 in Eq.©. The rea- 
son is found in a core radius which is substantially smaller 
than both the proton-proton distance as well as the dis- 
tance from their center of mass and the core, see table ED 
Combined with a smaller weight in Eq. Q on these terms 
for matter compared to charge radii this results in these 
larger charge radii. Recent measurements confirm that the 
charge radius is larger than the matter radius [571138] . If 
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Table 8. Root mean square (rms) radii in fm for the different 
bound states in 18 Ne with the three proton-core interactions. 
Charge and matter radii are shown in first and second row, re- 
spectively. They employ the root-mean square radius of 2.71 fm 
for the 16 O core. The only known experimental values for 18 Ne 
are the root mean square charge radius, 2.971±0.020 fm [37] . 
and matter radius, 2.81±0.14 fm [3H], for the ground state. In 
|37] we can also find FMD calculations that give 2.93 fm and 
2.70 fm for the charge and matter radii of the ground state. 
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we compare our values with FMD calculations [37] we can 
see that our matter radius is bigger and closer to experi- 
mental data. On the other hand the FMD charge radius 
is in better agreement with the experimental value. 

The average size can be distributed between distances 
of the different constituents, i.e. in the present case the 
proton-core and proton-proton distances. These root mean 
square radii for the computed bound states are given in 
table IH1 for the different potentials. These two-body dis- 
tances within the three-body system also follow the gen- 
eral trends of binding energy and angular momentum. The 
distance between the two protons is in all cases larger than 
the proton-core distance. This is because the proton-core 
attraction is decisive for the binding of all three-body 
states. For the weakly bound 3 + state this difference is 
substantially larger due to the small binding energy. The 
sizes show that the protons are located substantially out- 
side the surface of the core. This is obviously helping to 
decouple core and valence degrees of freedom, and validate 
the model assumptions in the treatment as a three-body 
system. In general the trends from the overall rms radii in 
table [5] are maintained. 

The average sizes in tables [8] and [9] are results of the 
probability distributions. They are shown in FigUfor the 
two + states as functions of the distances between the 
two protons and their center-of-mass and the core. Both 
distributions have a tail in the proton-proton distance ex- 
tending to about 10 fm. The fall-off in r CjPP seems to be 
faster reaching no more than about 5 fm. The ground state 
has two separated peaks around the points (r pp , r c ^ pp ) « 
(1.5,2.7), (5.0,0.8) (all in fm) where the last is much 
smaller and somewhat broader. In contrast the second + 
state has only one peak at (r pp ,r c ^ pp ) w (4.2,2.2). Since 
the angular momentum decompositions are rather similar 
these differences must arise from the interference between 
the adiabatic components. The prominent peak in 0^ is 



Table 9. For the different computed states in Ne, and the 
different proton-core potentials, root mean square distances 
(in fm) (rpp) 1 / 2 and (r C p) x ' 2 , where p and c denote an external 
proton and the core, respectively. 
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Fig. 4. Contour diagram for the probability distribution of 
0^ (upper) and 0^ (lower) states in 18 Ne. The square of the 
three-body wave function is integrated over the directions of 
the two Jacobi coordinates. 



in Oj moved to larger distances between the protons and 
the small peak is at the same time moved to somewhat 
larger distances in r c>pp . The result is that 0^" has one 
broad peak. 

The same pattern is found for the two 2 + -states with 
two peaks for 2^ at (r ppi r C)PP ) « (1.9, 2.7), (5.0, 1.0) and 
one broad peak for l\ at (r pp ,r CiPP ) w (5.0, 1.8). The po- 
sitions are also almost the same as for the + -states where 
the latter position tends to be at smaller distances. The 
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Table 10. Root mean square (rms) radii in fm for different 
negative parity resonances in 18 Ne with proton-core interac- 
tions III (for the 17/ state) or IV (for the hj, CP, 3j~, and 2~ 
states) . Charge and matter radii are shown in second and third 
column. We used the root-mean square radius, 2.71 fm, of the 
16 O ground state instead of the unknown value for the 3~ ex- 
citation. In the fourth and fifth column we show root mean 
square distances (in fm) (r pp ) 1//2 and (r cp ) 1//2 , where p and c 
denote an external proton and the core, respectively. 





(rl) 1 ' 2 


(r 2 ) 1 ' 2 


i r pp) ' 


(rc P ) 1/2 


ir 


3.00 


2.88 


5.80 


4.12 




2.86 


2.80 


4.36 


3.26 




2.77 


2.75 


4.11 


3.10 


sr 


2.72 


2.72 


3.64 


2.88 


2i 


3.01 


2.90 


4.37 


3.87 



4 + -state has one peak at (r pp ,r CtPP ) w (3.5,2.0) which is 
an almost equal sided triangle. The 3 + -state has one peak 
at (r pp ,r c . pp ) w (4.6,2.2). 

4.1.2 Unbound states 

The negative parity states are in principle all resonances 
but except for the lj~ 3^~ states they are all computed 
as bound states with respect to the 3~ core excitation. 
Therefore the radial moments are for these (3~ based) 
well defined and together with partial wave decomposi- 
tion characteristic for the structures. In table [10] we give 
root mean square radii of matter and charge together with 
distances between protons and core for these states. Here 
in order to calculate matter and charge radii we have used 
the same radii for the core. 

As for the positive parity bound states the charge radii 
are usually slightly larger than matter radii, see table [TU1 
In general the positive and negative parity states are al- 
most of the same size even though the binding to the 3~ 
core is larger by several MeV. Also the internal distances 
remain essentially the same with the proton distance as 
the largest. Not surprisingly the tendencies with binding 
energy and angular momentum follow the general rules 
explained for the positive parity states. Only the lj~ state 
of roughly zero energy has a sufficiently well defined radial 
structure to be included in table [TU] 

The probability distributions are all rather similar. 
Both l~-states have only one well-defined peak at (r pp , r CtPP ) 
w (1.9, 2.7), (3.5, 1.8). The first of these resemble the low- 
est + and 2 + state and the second more resembles the 
4 + state. We also find one peak for both 0~, 27/ and 37/- 
states at {r pp ,r C:PP ) w (3.1, 1.6), (1.9, 2.2), (2.3,1.6). For 
2^~ a second peak is indicated at the side of the main peak 
resulting in an intermediate structure between the two we 
have already seen for bound states (see Fig. [SJ. The 0~ 
state resemble the 4 + state with the two protons close to 
the core but rather far from each other. 

The probability distributions for the resonances built 
on the + core ground state are not well defined since their 
energies are above threshold and the radial wavefunctions 
therefore spread out to infinitely large distances. However, 




Fig. 5. Contour diagram for the probability distribution of 27/ 
state in 18 Ne. The square of the three-body wave function is 
integrated over the directions of the two Jacobi coordinates. 

they all have large amplitudes around the minimum in the 
corresponding adiabatic potentials. 

4.2 Electromagnetic decays 

The bound states can only decay electromagnetically. The 
corresponding observable transition probabilities are crit- 
ically depending on the structures. Thus they provide ex- 
perimental tests and we therefore compute the lifetimes for 
future comparison. The three-body states below 6.13 MeV 
built on the core-excited 3~ state are also bound states. 
They can therefore only decay by 7-emission to lower ly- 
ing states either by maintaining the same cluster structure 
or by -E3-decay of the core excited state. In both cases we 
can compute the electromagnetic transition probabilities. 
The selection rules determine the dominating transitions 
which can be of both electric and magnetic origin. 

The effective charge of the proton e e / / is for these es- 
timates determined by renormalizing the calculated E2- 
transition strength from first excited to ground state in 
17 F to the measured value of 25 Weisskopf units or 62.96e 2 fm 4 
|22j . This gives a value amazingly close to unity e e ///|e| = 
0.99 which is used in the three-body computations. 

The structures are sufficiently similar for the different 
interactions to allow estimates with only one potential for 
each state. We choose potential I and IV for positive and 
negative parity states, respectively. The selection rules de- 
termine the dominating transitions which can be of both 
electric and magnetic origin. Transitions within the same 
parity are then dominated by E2 or Ml emissions, and 
between different parity states predominantly by El. If 
El is forbidden the much weaker M2 or E3 transitions 
may determine the lifetime, but to this level of accuracy 
the neglected mixture of ground and excited core-states 
could contribute. 

4.2.1 Multipole operators 

The electric multiple operators arc defined as: 

A 

M^(EX) = e^Zi^Yx^fi) (10) 

i=l 
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where A is the number of constituents in the system, each 
of them with charge eZi, and where is the coordinate 
of each of them relative to the A-body center of mass. 
The electric multipole strength functions are defined 

as: 

B(EXJi^If) = ]T {(IfMflM^EX^Mi)] 2 

liM f 

= ^ n -K//||7W(^A)||/ i )| 2 (11) 

The M.^(EX) operator in Eq.(JTU]) can then be rewrit- 
ten as: 

A-2 

M^EX) = e Z,\r c + r'^Yx^r^+r',) 

2 

+ eY J Z 3 r*Y x ^{r j ) (12) 

i=i 

where the index i runs over the A — 2 constituents in the 
core, and j labels the two external nucleons. 

As in [5] we can rewrite into intrinsic core and external 
valence nucleon coordinates. The results turn out to have 
the form 

3 

M^EX) = e^Zi^Yx^ri) + M^{EX,c) (13) 

i=l 

A-1 k 

+ f\( k , m ^) r c Y k, m (f c )Mf,- m {E(X-k),c) 

k—1 m— — k 

where the two first terms refer to independent valence and 
core degrees of freedom, respectively. The last terms de- 
scribe simultaneous transitions of core and valence parti- 
cles where f(m, /i) is a well defined function of its indices. 
Then, since the only allowed core transition is A^ M (i?3), 
only the second term contributes to transitions between 
the two different core states. For transitions between the 
same core state the mixed terms may in principle con- 
tribute. However, these terms are accounted for by the ef- 
fective charge of the proton which was adjusted to describe 
the i?2-transition in 17 F. Thus also these terms should not 
be included. 

The magnetic multipole operator, A/f M (MA), has the 
opposite parity of M.^{EX) and give rise to much smaller 
rates for the same A. Thus M^(M2) only is active be- 
tween negative parity states where the E3 core transition 
is necessary in the present cases. The (M2) transitions are 
then forbidden. On the other hand (Ml) conserves par- 
ity, allows unchanged core structure, and may compete 
with (E2) transitions. We shall therefore only consider 
the (Ml) operator which is defined as: 

M » {M1) = 2ik{i YS^ s ^ + > ( 14 ) 

i 

where fi labels the spherical component of an operator, 
and the constants g s and gi depend on the constituent 



Table 11. i?2-transition strengths from first and second 2 + 
states to neighboring J + states expressed in units of e 2 fm 4 . 
The first two rows (upper part) contain measured values from 
[22] . The next two rows (central part) correspond to our calcu- 
lation. The last two rows (lower part) show computed values 
obtained in |20| . 



B 




o 2 + 


2+ 


4+ 


B(E2,2+ - 


*J+) 


48(5) 


2.7(13) 




43(6) 


B(£2,2+ - 


*J + ) 


1.8(9) 








B{E2,2+ - 


+ J+) 


21.53 


2.199 


2.445 


31.71 


B(E2,2+ - 


+ J+) 


0.322 


7.154 


2.445 


1.179 


B(E2, 2+ - 


+ J+) 


32.09 






59.26 


B(E2, 2+ - 


+ J+) 


0.084 









particles i. The magnetic multipole transition strength is 
defined as for the electric case (see Eq.(fTTj)'). 

(c) 

The core has a charge of 8 units suggesting that g\ = 
8. For the positive parity states where the core angular 
momentum is zero we use a vanishing effective spin g- 

(c) 

factor, i.e. <?J = 0. For the negative parity states where 
the core angular momentum is 3 the g-factor is unknown 
but also of little interest since the dominating decay prob- 
abilities are determined by other transitions. We use the 
free proton value of =5.586 and we use again an ef- 
fective proton charge gi p ^=0.99. The relevant transition 
operators are then defined and we can compute the ob- 
servable transition strengths. 



4.2.2 Transition strengths 

The S( J B2)-transition for 17 F from 1/2 to 5/2 amounts to 
about 63e 2 fm 4 corresponding to 25 W.u. and a width of 
1.6 x 10 -6 eV [22] . This is within 1% the same as computed 
from the relative two-body wavef unctions. This indicates 
that an effective charge of 1 should be used for the low- 
energy sequence of states in 18 Ne which all are dominated 
by E2 transitions. We collect in table QTJ the four possible 
B(E2) values from first and second 2 + to the other posi- 
tive low- lying parity states. As usual the reduced matrix 
element should be multiplied by 1/(2 E + 1) depending on 
initial and final states in the transition. The table values 
include this factor and reflect the chosen direction of the 
transition. 

For potential I the computed table values (central part 
of table ITT]) are all systematically smaller than the mea- 
sured results (upper part of table ITTj) , i.e. smaller by fac- 
tors of 2.2, 1.22, 1.36, respectively for the known 2^ tran- 
sitions in the third row of table [TT] For B(E2,2^ 0+) 
there are newer measurements resulting in 23 ± 4e 2 /m 4 
and 27±4e 2 /m 4 in |40) which are in better agreement with 
our values. However, later on the same author published 
in [UJ the results of a measurement of the life time of this 
excited 2+ state which is more consistent with the value 
we used from [22]. The computed 2\ — > Of transition is 
smaller by a factor 5.5 but on top of the varying experi- 
mental results the measured value in |22] is given with a 
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Table 12. E2 and Ml-transition strengths from the second 3 + 
state to lower-lying J + states expressed in units of e 2 fm 4 and 
e 2 /m 2 , respectively. The third row gives the ratio of partial 
widths for the two decay modes. The fourth and fifth rows 
shows the values from [20] 



B 


2+ 4+ 2+ 


B(E2,3+ -> J+) 
B(M1,3+ -> J + ) 
i^/An 


0.0596 24.0 4.27 
0.0240 1.76 x 10 -7 1.16 x 10~ 3 
1.35 x 10" 5 147 0.00250 


B(E2,3 + -> J+) 
B(Ml,3+ -> J+) 


0.107 5.89 14.4 
1.54 x 10" 3 3.84 x 10" 3 1.35 x 10" 4 



rather large uncertainty which could reduce this discrep- 
ancy to a factor of 2.7. The discrepancies are reduced if 
we correct all numbers by scaling the root mean square 
value of the ground state from the calculated 2.78 fm to 
the measured value of 2.81 fm as found in [35]. 

The remaining deviations are now within acceptable 
ranges for a model where core polarization is neglected. 
The largest discrepancy appears for the + state which 
might have the strongest influence from the lowest core 
excitation of the same quantum number, + . The phe- 
nomenological procedure to correct for that effect is to use 
an effective charge larger than unity which then accounts 
for influence beyond the single-particle degrees of freedom. 
Most of the transitions are substantially larger than cor- 
responding to one single particle unit. This usually is a 
signal of the need for an effective charge larger than unity 
which in turn implies that core degrees of freedom are 
important. On the other hand the transition probabili- 
ties are not consistent with collective vibrational motion 
as the spectrum of excited states otherwise could indi- 
cate. Furthermore, the main variation is picked up by the 
three-body model. 

With the radius of 1.2 A 1 / 3, fm the single particle Weis- 
skopf unit for B(E2) is 2.690e 2 fm 4 which indicates that 
the proton-core i?2-transition of about 25 W.u. is con- 
structed by substantially more than a simple single-particle 
transition. On the other hand this E2 transition in 17 F is 
reproduced with the proper wavefunctions and an effective 
proton charge of one. Then it is not unreasonable to ex- 
pect that the three-body system should be approximately 
describable without active core degrees of freedom as well 
or perhaps rather with appropriate effective charges from 
the two-body subsystem. 

The transition probability from 2\ to 2\ receives also 
a contribution from Ml. We find B(M1, 2+ 2+) = 9.7 x 
10~ 4 e 2 fm 2 which is smaller than the measured value of 
0.0017 ± 0.0007 e 2 fm 2 (0.088±0.038 W.u.) corresponding 
to a width of 9.5 x 10 -3 eV [22] . The small computed value 
has the inherent uncertainty arising from spin polarization 
which can lead to a substantial correction. In any case this 
decay seems to be dominated by Ml. 

The decay modes of lowest multipolarity for the 3 + 
state are E2 and Ml where the final state can be any 
of the three states shown in table [T2] The corresponding 
decay widths are given by A 2 = Ait /7b{E 1 /h,c) b B{E2) 



and A/i = 167r/9(i; 7 /^c) 3 S(Ml). The ratio of widths is 
then Aa/An = Q.03(E 1 /hc) 2 B(E2)/B(Ml). 

The E2 transition from 3 + is dominated by decay to 
the 4 + state. The E2 transitions to the two 2 + states are 
smaller by factors of about 6 and 30, respectively. These 
decay probabilities are essentially completely arising from 
single-particle proton transitions between orbits around 
the 17 F structure. The relative sizes correspond directly to 
the probabilities of the largest E2 allowed configurations 
in tables [13J [Hand [H 

The Ml transitions are also essentially due to proton 
transitions between orbits around the 17 F structure. They 
are all very small but still varying by orders of magnitude. 
The sizes strongly indicate that the dominating parts of 
these transitions are forbidden by Ml selection rules due 
to the single particle character of the operator. The deci- 
sive selection rules are Al x = A£ y = 0, and As x = As y = 
0. For the £-part also AS = and analogously AL — for 
the s-part. The 3 + state is from table [T51 seen to be dom- 
inated by (£ x ,£ y ,L, s x , S) = (2, 0, 2, 1/2, 1), (0, 2, 2, 1/2, 1) 
in the second and third Jacobi coordinate system. The se- 
lection rules then prohibits Ml transitions to the 4 + state 
between the dominating components, see table 1141 while 
much smaller components still contribute resulting in the 
small value in table [T2] The transitions to the 2 + states 
can only proceed via small components, see table [TBI and 
the contributions in table [T2] are consistently also very 
small. 

In any case the small Ml transition probabilities can- 
not be precise since admixtures of various kinds or changes 
of the already contributing components could change the 
numerical results completely. Assuming the values in ta- 
ble [T^] it is interesting to see that decays into 2+ and 2 J 
are dominated by Ml while the 4 + state preferentially 
would be reached by E2. 

In the Ml computations the effective values of the g- 
factors are rather uncertain and polarization effects could 
change them. The M 1-transitions are particularly sensi- 
tive to spin polarization which is able to change g s sig- 
nificantly. The orbital g-factors can be expected to have 
less uncertainty on a level similar to the effective charges. 
Unfortunately, no values are available from the two-body 
subsystem which otherwise could have provided the ba- 
sis for comparing two and three-body effective ^-factors. 
More precisely, in this way be able to give an answer to 
which degree of adjustment to two-body data is necessary 
to reproduce the three-body transition properties. 

The cluster model in [20] has a very different structure 
with only two-body components, either a+ 14 or p+ 17 F. 
A number of transitions between positive parity states are 
computed. We give the transition strengths in the cap- 
tion of table [IT] and in table UJ The B(E2,2+ -> Of) in 
[2"tj] is between our value and measurement, B(E2, 2+ — > 
4 + ) is larger than measured where our result is smaller, 
B(E2, 2f — > Of) in ref [2UJ deviate from experiment by a 
factor of 20 compared to our factor of about 5. The tran- 
sition B(M1,2% -> Of) is also further away from mea- 
surements than ours but both values are very small. The 
transition from the 3 + state is the focus of the two-body 



14 



J. A. Lay et al.: Three-body structure of low-lying 1 Ne states 



cluster model in [20 . No experimental values exist. For 
both E2 and Ml transitions to the 2 + states, [20] obtains 
larger values than in the present work and vice versa for 
the 4 + state. 

The negative parity states are more limited in their 
decay modes because they have to change three units of 
angular momentum on the core. This means that the tran- 
sition has to be of order E3 or higher. The only contri- 
bution is then seen from Eq. (fT5|) to be accompanied by 
the core-transition with B(E3) = 259.5e 2 fm 6 (13.5 W.u.) 
corresponding to a width of 2.60 x 10 -5 eV. This value 
has to be multiplied by the overlap of the valence wave- 
functions. Therefore the results are proportional to the 
overlaps discussed in subsection l3. 3.21 but in any case very 
small. Furthermore, any admixture of a + component 
in the dominating 3~ core structure of these resonances 
would determine the decay probabilities. 

These decays would proceed through the small admix- 
tures of core-state + in the negative parity states or by 
the presumably much smaller core-state 3~ admixture in 
the positive parity states. Let us assume a small mixing 
amplitude of e of core-state + . Then the decay could be 
of much smaller multipolarity and the rate therefore much 
larger. If El is allowed it would dominate but should be 
reduced by the amplitude e. Thus the E3 core transition 
from J~ has to compete with the El rate related to the 
matrix element: 

(0+(We); J+ = | J" ± l\\M(El)\0 + {core); J~). (15) 

The El rates should then be multiplied by e 2 indicating 
that even with a 0.1% admixture the i?l-transition would 
dominate. 



5 Summary and conclusions 

The properties of 18 Ne have been investigated assuming 
a three-body structure with an 16 core and two pro- 
tons. The aim is to establish these structures which are 
necessary ingredients in computations of momenta of pro- 
tons and 16 from two-proton decays. We use the hy- 
persphcrical adiabatic expansion method for the Faddeev 
equations. The proton-proton interaction reproduce low- 
energy scattering data. The proton-core two-body inter- 
action is adjusted to reproduce the bound s and d-states 
of the proton- 16 ( 17 F) subsystem. The negative parity 
two-body 17 F states are unbound but may still be used 
to constrain the p-wavc interaction. The p-waves arc then 
treated in different ways, i.e. first by using the same poten- 
tial as for s and d-waves and second by adjusting indepen- 
dently to the measured levels. Three different potentials 
are constructed with the ground state structure of the + 
structure for 16 O with different prescription to account for 
the Pauli blocking of the occupied core states. 

The unbound negative parity states are relatively high- 
lying. They appear to be more related to core excited 16 
states than to single-particle p- waves. We then use the 
16 3~ excited state as building block. The sequence and 
spacing of the low-lying 17 F negative parity resonances 



support this interpretation as essentially uncoupled struc- 
tures. The corresponding interactions are adjusted to re- 
produce the lowest-lying negative parity states of 17 F with 
an inert 3~ excited 16 0. We follow this conjecture and ex- 
ploit the consequence of uncoupled three-body structures 
built on these different 16 O structures. 

The first results are energies and structures of the five 
lowest-lying bound states of 18 Ne. The energies are repro- 
duced with the usual accuracy in such three-body com- 
putations of about 200 keV at least for our potential I 
and omitting the second + state (around 500 keV). We 
do not employ phenomenological three-body potentials to 
fine tune these energies. This is only necessary when pre- 
cise decay properties are requested. The second series of 
results are bound excited states built on the 3 _ core ex- 
cited 16 O state and decoupled from the first series of three- 
body states. The third series are higher- lying states arising 
from the + structure of 16 O but for potentials which only 
are adjusted to proton-core s and d-waves. The p-waves 
are present and contribute without the relatively strong 
attraction necessary to reproduce the p-resonance. These 
states have energies below a few MeV above threshold 
where several states are known in the measured spectrum. 
Furthermore in most cases these states are separated by 
several MeV from the states of the same angular momen- 
tum and parity arising from the 3 _ core excited 16 state. 
These structures may therefore still be essentially uncou- 
pled. 

The detailed structures of all these states are extracted 
and discussed in terms of the partial waves of the subsys- 
tems building these three-body states. These predictions 
would be ingredients in future tests involving transfer into 
these states or breakup or decay from them. We also pro- 
vide root-mean square radii for all states for both matter 
and charge and division into proton-proton and proton- 
core distances. Not surprisingly the distributions are all 
peaked at distances corresponding to the minima in the 
adiabatic potentials. In some cases also a smaller peak ap- 
pear nearby due to different contributing configurations. 

The transitions between these states are rather sen- 
sitive to their electromagnetic structures. For the bound 
states E2 transitions dominate and we compute them and 
compare if possible to measured values. In all cases we 
find very large decay rates although systematically lower 
than experiments by a factor between 1.5 and 2. However 
in the 2^~ — > 0^~ transition we have a much larger differ- 
ence with the experimental results although this value has 
a large uncertainty, as already mentioned, which could re- 
duce this difference. The Ml transitions between bound 
states only contribute between the two 2 + states, and 
although almost forbidden also between the 3 + and the 
2 + ,4 + states. The Ml transition rate between 2 + states 
is an order of magnitude smaller than the measured value. 
The widths for transitions between 2 + states and from 3 + 
to the 2\ and 2% states are dominated by Ml while 3 + 
to 4 + is dominated by E2. 

We compare to values from a two-body cluster model 
with very different structures of all the investigated states. 
For the lowest positive parity 0, 2, 4 states the _E2-transitions 
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are in both models within factors of about 2 from exper- 
iments. The transition from second 2 + to ground state is 
improved by a factor of 4 in our model, still missing about 
a factor of 5 compared to the measurement. Transitions 
involving the 3 + state are very different in the two models, 
and here especially the very small Ml-values. 

Electromagnetic Transitions between states related to 
different core structures are in the schematic model de- 
termined by the E3-transitions between the core states. 
However, in reality even very small admixtures of + core 
structure in the wavefunctions of dominating 3~ core struc- 
ture could easily enhance these transitions by many or- 
ders of magnitude by allowing lower multipolarity. Fur- 
thermore, such small admixtures would lead to three-body 
decay widths of these resonances completely dominating 
over the electromagnetic decays. These decays would then 
proceed through the small component of the + core struc- 
ture and decisive as determined only by the strong inter- 
action. 

In summary, we have in details investigated the three- 
body structure of the low-lying states of 18 Ne. The bound 
state energies are rather well reproduced by properties of 
the two-body subsystems. The resonances can not be com- 
pared directly to measured values due to lack of detailed 
information. These resonances fall in two groups related 
to two different structures (0 + and 3~) of the 16 core. 
The small experimental widths of many of the 18 Ne neg- 
ative parity resonances can be explained by their main 
structure as bound states with an 16 excited 3~ state. 
Their decay is then through small admixtures of the 16 
ground state in their wavefunctions. The E2-transitions 
between bound states are systematically smaller than the 
measured values but within the range of moderate values 
of effective charges and ^-factors. These investigations are 
necessary for computations of three-body decays of 18 Ne 
resonances but useful directly by providing detailed infor- 
mation about the structure of low-lying 18 Ne states. 

ACKNOWLEDGMENTS 

This work was partly supported by funds provided by the 
DGI of MEC (Spain) under Contract No. FIS2008-01301. 
J.A.L. acknowledges a Ph.D. grant from the University of 
Seville and C.R.R. acknowledges support by a predoctoral 
I3P grant from CSIC and the European Social Fund. 

APPENDIX 



A Wavefunctions Components 

We give the contributions larger than 0.1% of each Fad- 
deev component for a number of computed states as dis- 
cussed in section 13.31 Results for both the two different 
Jacobi coordinates are given for each state. 

References 

1. A.S. Jensen, K. Riisager, D.V. Fedorov and E. Garrido, 
Rev. Mod. Phys. 76 (2004) 215. 



Table 13. The same as table for the 2+ states in 18 Ne. 
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Table 14. 


The same as 


table [5] for the 4+ state in 18 r 
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100 


32.7 


33.4 


32.4 


4 





4 








60 


9.4 


9.0 


9.7 


4 





4 


1/2 





40 


0.1 


0.1 


0.1 


3 


1 


4 


1/2 





60 


0.9 


1.0 


0.9 


3 


1 


3 


1/2 


1 


60 


0.2 


0.2 


0.2 


1 


3 


4 


1/2 





40 


0.0 


0.1 


0.1 


1 


3 


3 


1/2 


1 


60 


0.1 


0.2 


0.2 


2 


2 


4 


1/2 





100 


31.8 


32.1 


32.6 


2 


2 


4 


1/2 


1 


60 


0.1 


0.1 


0.1 


2 


2 


3 


1/2 


1 


120 


66.7 


66.2 


65.9 
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Table 15. The same as table [5] for the 3 + state in 18 Ne. 



4 


p 

Cy 


L 


<5;e 


S 


-K-max 


Wi 


Wn 


Wm 


l 


1 


2 


1 


1 


120 


54.2 


54.1 


77.9 


l 


3 


2 


l 


1 


120 


19.4 


19.2 


8.2 


l 


3 


3 


l 


1 


40 


0.0 


0.1 


0.0 


3 


1 


2 


l 


1 


120 


19.8 


20.2 


8.2 


3 


1 


3 


l 


1 


40 


0.0 


0.1 


0.0 


3 


3 


2 


1 


1 


100 


6.6 


6.4 


5.7 


1 


1 


2 


1 


1 


100 


1.2 


1.1 


1.2 


1 


3 


2 


1/2 


1 


40 


1.8 


1.8 


1.4 


3 


1 


2 


1/2 


1 


40 


0.0 


0.1 


0.0 


2 





2 


1/2 


1 


120 


50.1 


50.0 


49.9 





2 


2 


1/2 


1 


120 


46.8 


40.9 


47.4 


2 


2 


2 


1/2 


1 


40 


0.1 


0.1 


0.1 


Table 16. 


The 


same 


as 


table [S] for the 0' 


" state in 18 Ne 


the 3" 


core 


excited state. 








4 


Cy 


L 


Sx 


S 


Kmax 


Wiv 






1 


1 


2 


l 


2 


100 


83.1 






1 


3 


2 


l 


2 


80 


7.7 






1 


3 


3 


l 


3 


60 


0.1 






3 


1 


2 


l 


2 


80 


7.7 






3 


1 


3 


l 


3 


60 


0.1 






1 


1 


2 


5/2 


2 


60 


0.1 






2 


2 


3 


5/2 


3 


80 


2.3 






2 





2 


5/2 


2 


100 


49.1 









2 


2 


5/2 


2 


100 


48.3 






1 


3 


2 


5/2 


2 


60 


0.1 






3 


1 


2 


5/2 


2 


60 


0.1 







Table 19. The same as table [5] for the 2 states in 18 Ne with 
the 3~ core excited state. 



4 


4 


L 




S 




Wiv 


1 


l 


2 


l 


2 


80 


19.8 
75.2 


l 


i 


2 


l 


3 


80 


9.6 
18.9 


l 


i 


2 


l 


4 


80 


2.6 
0.6 


2 





2 





3 


100 


24.8 
2.4 





2 


2 





3 


100 


43.2 
2.9 


2 





2 


5/2 


2 


100 


11.5 
34.8 


2 





2 


5/2 


3 


100 


25.8 
7.4 


2 





2 


7/2 


3 


100 


11.7 
7.9 


2 





2 


7/2 


4 


80 


1.5 
0.7 





2 


2 


5/2 


2 


100 


11.4 
39.0 





2 


2 


5/2 


3 


100 


7.6 

8.9 





2 


2 


7/2 


3 


100 


29.0 
4.7 





2 


2 


7/2 


4 


80 


1.5 
0.6 



Table IT. The same as table [5] for the first 3 X state in 18 Ne 
built on the 3~ core state. 



4 


Cy 


L 


Sx 


S 


K 


Wiv 
















3 


120 


97.4 




1 


1 





1 


3 


100 


1.5 




1 


1 


1 


1 


2 


60 


0.1 




2 


2 








3 


100 


1.0 













5/2 


3 


100 


42.0 













7/2 


3 


100 


55.7 




1 


1 





5/2 


3 


60 


1.1 




1 


1 





7/2 


3 


60 


1.2 




Table 


18. 


The 


same 


as 


table [5] for the 3; 


state 


the 4 


core 


ground state 


using potential I. 


4 




L 


Sx 


S 


K 


Wi 


Wm 





3 


3 








120 


54.8 


55.7 


1 


2 


2 


1 


1 


100 


31.0 


31.3 


1 


2 


3 


1 


1 


60 


1.0 


0.9 


2 


1 


3 








80 


11.1 


10.0 


3 





3 


1 


1 


60 


2.1 


2.1 


1 


2 


2 


1/2 


1 


80 


15.6 


15.7 


1 


2 


3 


1/2 





100 


27.9 


27.9 


1 


2 


3 


1/2 


1 


60 


2.2 


2.1 


2 


1 


2 


1/2 


1 


80 


17.4 


17.4 


2 


1 


3 


1/2 





100 


33.3 


33.1 


2 


1 


3 


1/2 


1 


60 


2.3 


2.2 





3 


3 


1/2 





60 


0.3 


0.3 


3 





3 


1/2 





60 


0.9 


1.1 
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